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Abstract 

Pion and kaon correlations in relativistic nuclear collisions are studied in the frame¬ 
work of boost-invariant, cylindrically symmetric hydrodynamics. It is investigated how 
the inverse widths, R ou t, R s \de, of the two-particle correlation functions in out- and 
side-direction depend on the average transverse momentum K± of the particle pair, 
the initial energy density eo, and the equation of state of the system. The QCD tran¬ 
sition leads to a time delay in the expansion of the system and consequently to an 
enhancement of the ratio R ou t/R s ide- This time-delay signal is found to be partic¬ 
ularly strong for large average transverse momenta K± ~ 1 GeV and initial energy 
densities accessible at RHIC, eo ~ 10 — 20 GeV fm -3 . Neutral kaon pair correlation 
functions, which are not influenced by final state Coulomb effects and less contami¬ 
nated by resonance decays than pion correlation functions, seem to be the ideal tool 
to detect this collective time-delay signature of the QCD transition. 


^Supported by DFG, BMBF, GSI, and DOE. 
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1 Introduction 


Lattice studies of quantum chromodynamics (QCD) exhibit a chiral symmetry restoring, 
deconhning transition around a temperature T c ~ 160 MeV (for (net) baryon-free matter) 
[lj. The only possibility to study this transition and the properties of the high-temperature 
phase of nuclear matter, the so-called quark-gluon plasma (QGP) under laboratory 
conditions is via relativistic heavy-ion collision experiments Q. 

Promising signals for the detection of the QGP emerge from the influence of the QCD 
equation of state (EoS) on the collective dynamical evolution of the many-particle system. 
Such signals have to be studied via relativistic hydrodynamics Q, since that is the only 
dynamical model which provides a direct link between collective motion and the EoS of 
strongly interacting matter. For instance, it was shown in 000 that the QCD transition 
softens the EoS in the transition region and therefore reduces the tendency of matter to 
expand. This has two important consequences. On one hand, the expansion of the system is 
delayed and its lifetime is considerably prolonged 1, 0 , 0 - On the other hand it was shown 
in H j9J that the reduced expansion tendency of matter leads to a reduction of the transverse 
directed flow in semi-peripheral heavy-ion collisions. 

In |l3 the hydrodynamic flow pattern was studied for spherically symmetric (so-called 
“fireball”) expansion as well as for cylindrically symmetric transverse expansion with lon¬ 
gitudinally boost-invariant initial conditions (so-called “Bjorken cylinder” expansion), both 
as a function of the initial energy density eo and a possible finite width AT of the QCD 
transition region. The QCD transition was shown to considerably prolong the lifetime of 
the expanding system as compared to the expansion of an (ultrarelativistic) ideal gas, pro¬ 
vided the initial energy density lies in a certain range of values. This range corresponds to 
AGS initial energy densities for the fireball geometry, and RHIC initial energy densities for 
the Bjorken cylinder geometry (for realistic values of the initial thermalization time scale 
t 0 ). Motivated by Pratt’s and Bertsch’s conjecture mm, it was shown in jl0|, that for 


both geometries the prolongation of the lifetime could be experimentally observable via the 
ratio -Rout/-Rside of the inverse widths of the two-particle correlation functions in out- and 
side-direction, since this ratio follows the behaviour of the lifetimes rather closely. 

In this paper we extend the investigations of |T0| for the Bjorken cylinder geometry in 
two ways. First, we study two-kaon in addition to two-pion interferometry. While pion 
correlation functions are contaminated to a large extent by resonance decays, one expects 
these to be of minor influence in the case of kaons |T3j . In addition, neutral kaon correlation 
functions are, in contrast to the pionic ones, not distorted by final state Coulomb effects 
which introduce an additional uncertainty in experimental data. 

Second, we now present a systematic study of -R ou t, -Rside, and their ratio as a function 
of initial energy density eo and of the mean transverse pair momentum K± at midrapidity. 
As in IdJ, we take an EoS with a first order phase transition (AT = 0) as well as an EoS 
with a finite width AT = 0.1 T c (as an upper limit deduced from QCD lattice calculations 
[IJ), and compare the enhancement of -R ou t/-Rside relative to the ideal gas case. 

We find that the excitation function of -R ou t/-Rside for kaons again mirrors closely that of 
the lifetime of the system, and has a maximum around the same initial energy density eo 
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as for pions. Due to kinematical reasons the ratio is, however, somewhat smaller for kaons 
with small average transverse pair momentum K± < ttik = 497 MeV as compared to the 
pion case. On the other hand, for high K± ~ 1 GeV, it is nearly of the same order as that 
for pions. 

We also study separately the dependence of R out and R S id e on e 0 and Ii±. As expected, 
these quantities in general decrease with increasing K at fixed e$. An important exception 
is the situation around values of eo ~ 10 — 20 GeV fm -3 , where the time-delay signal is most 
pronounced. For these initial energy densities and in the case of a first order transition, 
AT = 0, we find that R 0ut actually increases as a function of K±, both for pions and kaons. 


We analyze this behaviour in detail studying the emission function and its variances |14 


The outline of the paper is as follows. In Section 2 we briefly discuss the hydrodynamic 
model, the EoS used, and the basic formulae to derive single inclusive particle spectra and 
two-particle correlation functions from the hydrodynamic solution. In Section 3 we present 
our results. A summary and conclusions are given in Section 4. An appendix contains 
technical details on the calculation of the variances of the emission function. Natural units 
U = c = k B = 1 are used throughout this work. 


2 Hydrodynamics and Particle Spectra 

2.1 The hydrodynamic model 

Hydrodynamics is equivalent to local conservation of energy and momentum, 

d„T ,M ' = 0 , 


( 1 ) 


plus additional conservation equations for the (net) 4-current of any conserved charge in 
the system. The latter will not be considered here, since our focus of interest, the central 
(rapidity) region of ultrarelativistic heavy-ion collisions, is expected to be essentially free of 
(net) charges due to the limited nuclear stopping power. 

In order to solve the relativistic hydrodynamic equations (T) we make the assumption of 
local thermodynamic equilibrium, the so-called “ideal fluid approximation”. This approxi¬ 
mation, although commonly employed throughout the literature, is rather restrictive, since 
the fluid is certainly out of thermodynamic equilibrium in the early and late stages of a 
heavy-ion collision. On the other hand, a commonly accepted, causal theory of dissipative 
relativistic hydrodynamics, which would account for non-equilibrium phenomena, does not 
yet exist 0. For the moment we have to rely on the assumption that (local) thermaliza- 
tion does indeed occur after the initial pre-equilibrium stage which is dominated by (hard) 
parton-parton scattering processes (for collisions at RHIC energies). Arguments in favour 
of this assumption are made in (T5j . We then describe the system’s evolution from this point 
on via ideal fluid-dynamics, until dissipative effects in the final stages prior to decoupling 
become large and the assumption of local thermodynamical equilibrium again breaks down. 

In the ideal fluid approximation, the energy-momentum tensor reads []H| 


= {e+p) u 


pg 


fils 


( 2 ) 
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Here, e, p are energy density and pressure in the rest frame of the fluid, while = 7 (1, (3) 
and g^ v = diag(+, —, —, —) are the 4-velocity of the fluid and the metric tensor, respectively 
((3 is the 3-velocity, 7 = (1 — j3 2 )~ 1 / 2 is the gamma factor). 

Our model for the dynamics of central ultrarelativistic heavy-ion collisions is the so- 
called “Bjorken cylinder” expansion [17], i.e., it assumes cylindrical symmetry perpendicular 
to the beam axis (taken to be the z—axis throughout this work) and a scaling solution for the 
longitudinal expansion ( j3 z = z/t). Note that due to the details of the mini-jet production 
mechanism, these assumptions about the symmetries of the system are not valid per single 
event |T 8 |. On the other hand, they seem quite reasonable on an event-average basis. 

A consequence of longitudinal scaling is boost-invariance along the beam axis. It is 
therefore sufficient to solve the system of four hydrodynamic equations ([!]) at z = 0 , where 
it reduces to 


d t E + d r [(E+p)(3] 
d t M + d r ( M(3 + p) 


^ + i 

r t , 


'P + 1 ' 

r t 


(E + p) , 
M . 


(3) 

(4) 


Here, (3 is the radial component of the transverse velocity, and we used the definitions 
E = T 00 , M = T 0r (r as index for the radial component of the energy-momentum tensor). 
In order to solve the (effectively (l+l)-dimensional) system ([|[4|), we apply the same scheme 
used in ITO 


i.e 


., we employ Sod’s operator splitting method [BJ to account for the source 
terms on the right-hand side, and the numerically well-tested relativistic HLLE algorithm 
H |2(J for the one-dimensional hydrodynamic transport on the left-hand side of d®. The 
HLLE scheme was shown to yield excellent results in particular for equations of state with 
phase transitions [Hj. pTJ. 


2.2 The equation of state 

As mentioned in the introduction, the advantage of hydrodynamics to model the dynamics 
of relativistic heavy-ion collisions is that this is the only approach which directly relates 
observable collective flow phenomena to the EoS of hot and dense nuclear matter. Formally, 
specifying an EoS closes the coupled system of hydrodynamic equations. In the particular 
case of vanishing conserved (net) charges considered here, the system (||.[f[) is closed by an 
EoS of the form p(e). In the ideal fluid approximation the emerging flow pattern is completely 
determined by the choice of this EoS and the initial conditions. 

I 11 the present work we study hydrodynamic expansion for three different equations of 
state. The simplest one is the EoS of an ultrarelativistic ideal gas, 

P(e) = c 2 s e , (5) 

where c s = l/y/3 is the velocity of sound. This is commonly taken to be a reasonable first 
approximation for a hadronic matter EoS in the central region of an ultrarelativistic heavy- 
ion collision, since here the system is hot, T > m^, and thus dominated by pions, which are 
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the lightest hadrons. The EoS ([]) has a constant pressure gradient dp/de = c 2 = 1/3 for all 
energy densities. 

To study the effects of a (phase) transition to QGP on the dynamics of the system, we use 
the following parametrization of the entropy density as function of temperature 00 . 01 . 


s 

Sc 


(T) 


T 

Tr 


\ + tanh 


cIq + d 


H 


T —T r 


AT 


( 6 ) 


Here, s c = (dq + du)^ 2 T ^/45 is the entropy density at the critical temperature T c , and 
dq and g?h are the number of (massless) degrees of freedom in the QGP and the hadronic 
phase, respectively (we fix dq = 37 and du = 3 throughout this work, for the effect of 
other choices on the dynamics of the system, see |]lT3|] ). The function p(e) then follows from 
thermodynamical relations and is used in the hydrodynamic model in tabular form. Eq. 
([]) accounts for a possible finite width AT of the transition. Present QCD lattice data [JT|J 
indicate 0 < AT <0.1 T c . 

In the limit AT = 0, eq. ((6) exhibits a first order phase transition between an ultrarcl- 
ativistic gas of hadrons described by the EoS (||), and a (net) baryon-free QGP consisting 
of gluons and u and d quarks (and antiquarks) described by the well-known MIT bag EoS 
Energy density and pressure in the QGP phase are then e q = dqTT 2 T A /30 + B, p q = 
dB)/3, while the corresponding quantities in the hadronic phase are 6 h = d^^T 4 /30, 
Ph = eh/3. B = \ [(d Q — du)/(dq + g?h)] T c s c is the MIT bag constant, if one measures en¬ 
ergy densities in units of the critical enthalpy density w c = e c + p c = T c s c (p c is the pressure 
at T c ; for our choice dq/du = 37/3 and T c = 160 MeV, T c s c ~ 0.75GeVfm -3 in physical 
units). Besides this first order transition with AT = 0, we also consider a transition with 
a finite width AT = 0.1 T c to cover present uncertainties in the QCD EoS. Note, however, 
that a finite width would also naturally emerge in any hnite system, even if the transition is 
of first order in the thermodynamic limit [|24] . 

The three different equations of state are displayed in Fig. 1. Part (a) shows the entropy 
density divided by T 3 and (b) the energy density divided by T 4 as functions of temperature. 
For the equations of state with a transition, at T = T c both quantities exhibit the strong 
increase (or even a discontinuity) observed in lattice calculations. Part (c) shows the pressure 
and (d) the velocity of sound squared, c 2 s = dp/de, as functions of energy density. While for 
the ideal gas EoS (dashed lines) c 2 = 1/3 is constant for all energy densities, this is different 
for the EoS with a first order phase transition, AT = 0 (full lines), and the EoS with a 
transition of hnite width AT = 0.1 T c (dotted lines). The pressure gradient c 2 s = dp/de 
vanishes for AT = 0 as long as matter is in the mixed phase, en = eh(T c ) < e < e q (T c ) = 6 q. 
For parts of the system, which pass through the corresponding region of energy densities, 
the tendency to expand is then considerably reduced compared to the ideal gas case, which 
results in a prolonged lifetime of the system B, 0 0 0 . For the transition with hnite 
width AT = 0.1 T c , the pressure gradients are hnite, but nevertheless smaller than for the 
ideal gas. Thus, again the system’s lifetime will be prolonged B 0 . 
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2.3 Particle spectra and interferometry 


In a heavy-ion collision, the ideal fluid approximation breaks down in parts of the system 
which become so dilute that the collision rates between the particles are smaller than the 
(local) expansion rate of the system and thus local thermodynamic equilibrium can no longer 
be maintained. These parts of the system decouple (“freeze out”) from the hydrodynamic 
evolution and, eventually, the corresponding particles freely stream to the detectors. In a 
rigorous sense, on one hand they can no longer be described in the hydrodynamical frame¬ 
work, while on the other, the energy and momentum carried away by these particles might 
influence the hydrodynamic evolution of the rest of the system. 

Although attempts have been made |25|, [26| , up to now no commonly accepted, consistent 
description of this so-called “freeze-out” process exists for the hydrodynamical description of 
heavy-ion collisions. The common way to treat this problem is to solve the hydrodynamical 
equations in the whole forward light cone and a posteriori identify those parts of the system as 
decoupled, which have cooled below a certain freeze-out temperature Tf or whose density has 
dropped below a freeze-out density rif. Final particle spectra are then calculated along the 
freeze-out hypersurface, £(:r), defined by the condition T(x) = Tf, or n(x) = rif, respectively. 

Here we also follow this prescription and employ the approach of Cooper and Frye 
to calculate the particle spectra (see e.g. [|28| for a (1 +1 )-dimensional application). In our 
case of vanishing (net) charges, the particle density is a function of temperature only, and 
thus freeze-out at constant density is equivalent to freeze-out at constant temperature. The 
freeze-out temperature is here taken as a free parameter in the range (0.7—1) T c ~ (110 — 160) 
MeV. 

The single inclusive momentum distribution of “frozen-out” particles, i.e., particles whose 
world lines cross the 3-dimensional hypersurface £ in 4-dimensional space-time is then given 
by 0 


dlV 


d 


d 3 k (27r) : 


■JdE-kf(k- u/T ) 


(7) 


In our case of pions or kaons, f(z) = (e 2 — l)^ 1 is the Bose-Einstein distribution function, 
dis the normal vector on the (local) hypersurface element £ M , and the integration runs 
over the complete freeze-out hypersurface £. Note that, while T is per our definition of 
freeze-out constant along £, the fluid velocity u M is not. 

The Cooper-Frye formula (0) is strictly valid only along space-like parts of the hyper¬ 
surface £, i.e., where the normal vector dE^ is time-like. For time-like parts of £ (space-like 
d£ M ), however, d£ • k can in principle become negative, corresponding to particles which 
re-enter the fluid and thus do not decouple from the system. From the physical point of 
view, this “negative” number of particles should not enter the single inclusive spectrum of 
frozen-out particles on the left-hand side of (^). Therefore, modifications of the Cooper- 
Frye formula have been proposed [^5|, ^9], which correct for this. They shall, however, not 
be employed here, because for the hydrodynamic solutions studied in the following we found 
the collective flow across the freeze-out surface to be rather strong. The maximum of the 
single-particle distribution function / is therefore shifted in a way which renders the number 
of particles with given 4-momentum that actually re-enter the fluid negligible. We shall 
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quantify this statement below (cf. discussion of Figs. 7, 8). 

The two-particle correlation function of two identical particles with momenta ki, k 2 is 
defined as the ratio of the two-particle probability P(ki,k 2 ) to the product of the single¬ 
particle probabilities for uncorrelated particles, 


C*2 (k x , k 2 ) = 


P( k ls k 2 ) 
P(ki)P(k 2 ) 


( 8 ) 


We follow the method outlined in MM M to compute (|D • Introducing the average 4- 
momentum of the particle pair = {ki-\-k 2 ) / 2, and the relative 4-momentum = k d — k 2 , 
and assuming that the particle source is chaotic and sufficiently large, the two-particle 

113 - 0 


correlation function for bosons can be written as 


C 2 (k|, k 2 ) — 1 + 


(4/ s d V-K exp [iq-Y] f (K ■ u/T) 


E 1 [dA^/d 3 ki] E 2 [dAT/d 3 k 2 ] 


(9) 


For the Bjorken cylinder geometry, the number of independent variables can be reduced 

m 


due to rotational symmetry around the beam axis. As in 


we focus on the correlation 
We choose the coordinate 


functions of particles emitted at midrapidity (K z = q z = 0) 
system such that K = (PA, 0,0), q ou t = (Qout,0,0) for the component of q parallel to K, 
and q si de = (0, g si d e , 0) for the component of q perpendicular to K. With these assumptions, 
C 2 (/l, g ou t, <?side) becomes a function of three independent variables only. 

The so-called out-and side-correlation functions are now defined as 


C*2,out (//out) — C ‘2 (Pj_, Qout; 0), 

Cgsick^Qside) — ^2 (Pj_, 0, (/side) • 


( 10 ) 

(ii) 


We then define the inverse width of the correlation functions in out- and side-direction by 


Pont — 


Qout 


Pside — 


Qside 


( 12 ) 


where q* ut is determined by C 2!OUt (g* u t) = 1.5, and q* ide analogously. The inverse widths ( |T2| ) 
are a qualitative measure for the duration of particle emission and the transverse size of the 
source. Their quantitative values, however, depend on the initial size of the system, Rq, on 
the transverse momentum K±, as well as on the details of the hydrodynamic flow pattern 
at freeze-out [II], O, [H |. As argued in [ID], these dependences should in principle largely 
cancel out if one considers the ratio of the inverse widths, P out /P S id e - In 13 it was also 


shown that, for the spherical expansion as well as for the Bjorken cylinder expansion, this 
ratio is a good qualitative measure for the lifetime of the system. The explicit form of (]9]) 
for the Bjorken cylinder geometry was given in Appendix B of |l(| and shall not be repeated 
here. 

For the side-correlation function, the single-particle spectra in the denominator of (|9|) 
are equal (because of rotational symmetry, the single-particle spectrum can only depend 
on the modulus of the transverse momentum, which is equal for the two particles, k 2 = 
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(K i, q s Me/2, 0), k 2 = (K±, —g si d e /2, 0) for pairs at midrapidity), while for the out-correlation 
function, they are approximately equal as long as the modulus of the relative momentum, 
g out , of the particle pair is small compared to K± (note that ki = (AT + g out /2,0,0), k 2 = 
(AT — g out /2, 0, 0) for pairs at midrapidity). Furthermore, in this case K° = (h® + )/2 ~ 

\J.K 1 + m? = Ex, and the correlation function can be approximated as ||T4| 


C 2 (K,q)~l + 


/ d 4 x exp [iq ■ x] S[x, K ) 


/d 4 a:A(a:,A') 


= 1 + | (exp [iq ■ x] 


where 


S(x, K) = -A f dE • K f(K • u/T) S ft \x - S) 
[An) 6 J s 


(13) 


(14) 


is the so-called emission function of the source, and the right-hand side of (|T3|) defines the 
average (•). In the following we shall make the Boltzmann approximation for the single¬ 
particle distribution function, f[z ) ~ e~ z , which usually holds to good accuracy, especially 
for large average pair momenta AT- 

In the case that the emission function is Gaussian in space-time, the correlation function 
can be written in the form (again, we focus on midrapidity pairs, K z = q z = 0) HIT 


C' 2 (AT, gout, <?side) ~ 1 + exp 


-ft 2 a 2 

1 ’'side uside 


p 2 „2 

-^out y out 


(15) 


where the ( Kj_ -dependent) so-called “side-” and “outwards radii” Aside, A ou t of the correla¬ 


tion function are related to the variances of the emission function via 14 


A; 


2 _ 


L out 

7+2 

-^side 


= <^ 2 ) 
= (?) 


2 Px(xt) + (3 2 k (t 2 ) 


(16) 

(17) 


Here (3x = K^/Ex is the (approximate) transverse velocity of the particle pair and x = 
x — {x) ( t , y are defined analogously). Note that, in the case of a Gaussian source leading 
to the correlation function (|15D, the inverse widths defined in (|l 2 l) are related to the radii 
(|T^,|T^) via A 0 ut,side — 1.2 A 0 ut,side- As will be seen in the next section, this formula holds 
to astonishingly good approximation in the cases considered here, even when the emission 
functions do not resemble Gaussians. The eqs. SB therefore allow for a reasonable qua¬ 
litative and quantitative understanding of the shape of the source or the emission function, 
respectively. The explicit formulae for the variances required to calculate A 0 ut,side f° r the 
Bjorken cylinder expansion are given in the appendix. 

We finally note that other parametrizations of the correlation function, such as the Yano- 
Koonin-Podgoretskii approach |32|], may be more advantageous than the standard (Pratt- 
Bertsch) parametrization from the point of view of obtaining a better quantitative estimate 
of the variances of the emission function [T|]. For studying the qualitative behaviour of 
the proposed time-delay signature, however, the standard parametrization proves to be 
sufficient. 


























3 Results 


For the Bjorken cylinder expansion one has to specify an initial (proper) time To for the 
expansion in order to avoid the divergence at t = 0 in eqs. (03). Physically, this time 
corresponds to the thermalization time scale prior to which initial, non-equilibrium parton- 
parton scattering processes dominate the dynamical evolution. Only for times r > To, 
dynamics can be described by means of (ideal) hydrodynamics. According to the “hot-glue” 
scenario 0 , r 0 ~ 0.5 frn for central Au+Au-collisions at RHIC. Since the transverse radius 
of the hot central zone in such collisions is of the order R 0 = 5 frn, we fix To = 0.1 i? 0 in the 
following. The dependence of the hydrodynamic flow pattern on varying To as well as the 
initial energy density e 0 and the EoS was discussed in detail in |Rj (cf. Figs. 7-9 therein) 
and shall not be repeated here. 

In particular, it was found in |T(|] that, for fixed average pion pair momentum K^_ = 300 
MeV and for an initial energy density eo = 18.75 T c s c , in the case of the QCD transition 
the experimentally measurable ratio R ou t/R S ide was considerably enhanced over the ideal gas 
case, reflecting the delay in the expansion of the system due to the softening of the EoS in 
the transition region. This is once more shown in Fig. 2 (thin lines) for the hydrodynamic 
expansion with a first order phase transition (AT = 0, dq/dn = 37/3, m n = 138 MeV). 
The three curves correspond to different freeze-out temperatures Tj = 0.7 T c (full lines), 
Tf = 0.9 T c (dotted lines), and Tf = T c (dashed lines). In this work we study this time-delay 
signature of the QCD transition not only as a function of eo, but also its dependence on K±, 
both for pions and kaons and for the three different equations of state discussed in Subsection 
2 . 2 . 

Note that all correlation functions in this work are calculated taking the finite mass of 
the respective particles into account, but for the fluid evolution, hadron matter is considered 
to be an ultrarelativistic ideal gas. While this appears inconsistent from a rigorous point of 
view, quantitatively we do not expect a finite pion mass in the EoS to severely affect the 
hydrodynamical evolution prior to decoupling, since the freeze-out temperatures considered 
here are around or even above m n . Kaons in the hadronic phase are, on the other hand, 
suppressed by their larger mass, mx > T c . Therefore, the effect of these particles on the 
hydrodynamic evolution is expected to be small. This view is supported by the results 
of [[33]. Moreover, small differences in the hydrodynamic evolution tend to have an even 
smaller impact on the correlation functions (which represent in a certain sense averages over 
the dynamical evolution) and thus on R out and i? S id e - 

The thick lines in Fig. 2 represent our new results for kaons (uik = 497 MeV, K± = 
300 MeV). For these heavier mesons we also observe a local maximum in the excitation 
function, which again mirrors closely the prolongation of the lifetime of the system around 
eo — 18.75 T c s c . Such a behaviour is not seen in the ideal gas case, where i? ou t/-R S ide varies 
rather smoothly (for the sake of clarity, we do not show results for the ideal gas in Fig. 2; for 
that case see, for instance, Figs. 3, 5 below and Fig. 17 (a) in Ref. [|T0|] ). Finally, note that 
for K l = 300 MeV, the enhancement is nearly independent of the freeze-out temperature 
both for pions as well as for kaons. 

For K l = 300 MeV the absolute magnitude of i? 0 ut/-R S ide is about a factor of 2 smaller for 
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kaons than for pions. The reason is the higher mass of the kaons compared to pions which 
leads, for the comparatively small value of K± = 300 MeV, to a considerable difference in 
the velocity of the kaon pair, /3k — 0.52, as compared to a pion pair with the same K±, 
/3k — 0.91. In Figs. 7, 8 below we shall explicitly show that the radii (|T6| , |T7| ) are reasonably 
good approximations for the actual inverse widths of the correlation functions. The most 
important term in ([T73| ) turns out to be (3\(i 2 ). Since /3\ ~ 0.27 for the kaon pair, but 
(T 2 K ~ 0.83 for the pion pair, even if the time variances ( t 2 ) of the emission functions for kaon 
and pion pairs were roughly equal, the outward radius would turn out to be much smaller 
in the kaon case. 

In fact, the time variance for the kaon emission function is in general smaller than that 
of the pion source. In order to understand this, note that the maximum of the emission 
function is typically located near points of the freeze-out hypersurface where the particle pair 
velocity /3k matches the fluid velocity (3 , since this minimizes the argument of the exponential 
term exp [—K ■ u/T] in the emission function. For a larger particle mass, however, the 
emission function decreases faster around this maximum, which in turn reduces its variance 
as compared to the case of a smaller mass particle. Since this effect, however, also influences 
the spatial variance of the emission function, we expect it to (roughly) cancel out in the ratio 
-Rout/-R S ide- Thus, the suppression of i? ou t/-R S idc f° r kaons as observed in Fig. 2 is mainly due 
to the smaller kaon pair velocity. 

On the other hand, following this argument we expect that for increasing values of Ii±, 
for which the pair velocity for kaons is closer to the causal limit (and to that of pions), 
i? ou t/i?side should approach similar values as for pions. This is indeed seen to be the case in 
Fig. 3 where we plot R 0 ut/R S idc f° r K± = 700 MeV and AT = 0. As in the pion case, there 
is a strong enhancement of about a factor of 2 over the corresponding ideal gas case (thin 
lines) for all three freeze-out temperatures. On the other hand, the dependence of R out /R s ide 
on the freeze-out temperature is larger than at K± = 300 MeV (cf. Fig. 2). The time-delay 
signal is strongest for late freeze-out (Tf ~ 0.7 T c ). 

To study the influence of K± on the correlation signal in more detail, we present in Fig. 
4 the ratio i? ou t/-Rside as a function of K± at fixed eo = 18.75 T c s c , i.e., where the time- 
delay signal is maximized, and for the case of a first order transition in the EoS. For pions 
(thin lines) as well as for kaons (thick lines) we observe a strong increase of the ratio with 
K l- While for pions, R 0 ut/R S ide is of order 3 already at low momenta, the ratio for kaons 
increases only gradually from around 1 at K± = 100 MeV, on account of the above discussed 
suppression due to the smaller kaon pair velocity. As expected, for high momenta, K± ~ 1 
GeV, R out /Rside for kaons approaches values comparable to those for pions, since both pion 
and kaon pair velocities approach unity. 

It is also seen from Fig. 4 that the observation made in [10] that the time-delay signal 
does not depend strongly on the freeze-out temperature is not valid in general and was 
only due to the specific choice K± = 300 MeV in that work. The sensitivity of the value of 
Rout/Rside on the freeze-out temperature increases with K±, such that R ou t/R S ide is the larger 
the smaller the freeze-out temperature is and may vary up to 20% for freeze-out between 
Tf = 0.7 T c and Tf = T c . 

In Fig. 5 we present the main result of the present work: shown are surface plots of 
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-Rout/-Rside in the (eo — Jij_)-plane. The freeze-out temperature was fixed at Tf = 0.7 T c . Fig. 
5 (a) shows the ratio of inverse widths for pions and AT = 0 (solid lines) versus the ideal gas 
case (dotted lines). In part (b) the same is shown for AT = 0.1 T c , while (c) and (d) show 
the corresponding results for kaons. As one observes, the overall behaviour of R ou t/Rside is 
rather similar for pions, Figs. 5 (a,b), as for kaons, parts (c,d). In accordance with Fig. 4 and 
for the (kinematic) reasons discussed above, the ratio at low K± is smaller for kaons than 
for pions. Also, the enhancement over the ideal gas case is less pronounced in the case of a 
smooth transition, AT = 0.1 T c , than for a first order phase transition, for reasons discussed 
in detail in Ref. [|(| 


The most important property of Fig. 5 is to allow to identify the range of initial energy 
densities and transverse pair momenta required to maximize the time-delay signature. For 
our choice of r 0 = 0.1 R 0 ~ 0.5 fin, the signal is maximized for initial energy densities around 
e 0 ~ 20T c s c ~ 15GeVfm -3 . For smaller (larger) r 0 , the respective values for e 0 increase 

10]. These values are well above the soft region of the EoS, 
|TU|, this is due to the strong longitudinal 


(decrease), as explained in 
cf. Fig. 1 (d). As discussed in detail in Ref. 


motion in scaling hydrodynamics and opens the possibility to observe the time-delay effect 
at the RHIC collider. Moreover, Fig. 5 confirms that, for sufficiently high K±, the signal is 
as strong for kaons as for pions. Kaon interferometry is, however, advantageous from the 
experimental point of view, since kaon yields in general are less contaminated by resonance 
decays, while neutral kaon correlations in particular are not subject to Coulomb distortions 
from final state interactions. 

In the following, we want to elucidate the behaviour of R ou t/R S ide in greater detail by 
studying the inverse widths separately as functions of e 0 and K±. Fig. 6 shows (a,c) R out 
and (b,d) R si d e both for (a,b) pions and (c,d) kaons for the case of a first order phase 
transition, AT = 0. As one observes, for fixed K±, R si d e increases slowly with e 0 . This 
behaviour is naturally explained by the fact that higher initial energy densities drive the 
transverse expansion of the system more strongly (cf. Ref. [10] for the explicit hydrodynamic 
solutions), which increases the transverse size of the source. On the other hand, for fixed eo, 
.Rside slowly decreases with increasing K±. To understand this quantitatively, recall that the 
emission function ([14]) has a maximum where K ■ u is minimized. For larger pair momenta 
K l (larger pair velocities Pk) the emission function S(x,K ) ~ exp [—K ■ u/T] will drop 
faster in the region around the maximum, which in turn makes the source appear smaller. 

On the other hand, the eo-dependence of R ou t shows again clearly the time-delay sig¬ 
nature of the QCD transition. An unexpected feature is, however, that R out increases with 
K± at fixed eo, at least for the case AT = 0 and for initial energy densities below about 
100 T c s c . According to the above argument, R out is expected to decrease as well, and the 
time-delay signature of Fig. 5 solely due to the fact that the decrease of R out is slower than 
that of -R S ide- 

In order to understand why this is not the case for AT = 0, let us first consider a situation 
where this expectation is actually fulfilled. Let us consider the hydrodynamic expansion in 
the case of a smooth transition AT = 0.1 T c and for an initial energy density eo = 18.75 T c s c 
(cf. Fig. 9 (c) of Ref. ]TD[). In Fig. 7 (a) we show the corresponding isotherms in the (t — r)- 
plane (this figure is identical with Fig. 9 (d) of [[b]]), part (b) shows the flow velocity P along 
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the freeze-out hypersurface, i.e., in our case the isotherm with T = 0.7 T c in Fig. 7 (a). Since 
the mapping of the isotherm on time is unique we prefer to plot (3 as function of t rather 
than as function of r, which, considering the shape of the isotherm, is obviously not unique. 

In Fig. 7 (c) we show the emission functions S(x,K) for pions, integrated over (j) and rj 
(i.e., more precisely, the integrand of eq. (|29|) in the appendix), as a function of time along 
the freeze-out isotherm for three different values of the average transverse pair momentum, 
K l = 300, 700, and 1000 MeV. As one expects, for larger pair momenta/velocities the 
maximum of the emission function shifts to regions of higher transverse velocity on the 
hypersurface, i.e., to earlier times. Also, in confirmation of the above argument, the emission 
function becomes narrower, which reflects in a decreasing R% ut , as shown in Fig. 7 (e). As 
expected, Rg ide decreases with K±, too, although this cannot be seen directly from the 
presentation of the emission function in Fig. 7 (c). 

Further understanding of this result can be gained by studying the radii (p~6|jr7|) . Accord¬ 
ing to eq. (0). consists of three different terms. The dashed line in Fig. 7 (e) shows the 
actual time variance ( t 2 ) of the emission function. In accordance with the above arguments, 
and in agreement with Fig. 7 (c), this quantity decreases with Ii±. For the dotted line, 
the time variance was multiplied by For the light pions, (3\ is close to unity already 
for moderate K± ~ 300 MeV and the difference between these two curves rapidly vanishes 
with increasing K±. For the long-dashed line the term —2 (3k {xt) was added. Obviously, 
the x — f-correlation of the emission function is negative and about 50% of the value of 
Px{t 2 ). The dash-dotted line is the full result for Rl ut . As one notices, the variance in 
^-direction is only a minor correction to the other two terms. It is astonishing how well 
the approximate relations ([l(|jT7]) reflect the qualitative and quantitative behaviour of the 
inverse widths of the correlation functions, although the shape of the emission function is 
not Gaussian. (Note that all variances appearing in eqs. ([L6|,|T7]) have been multiplied with 
a factor 1.44, to account for the difference in the definition of inverse widths and radii, as 
discussed at the end of Section 3.) 

In Figs. 7 (d,f) we show the corresponding results for kaons. One observes that the 
emission function in the case K± = 300 MeV is now somewhat narrower in time direction 
than in the pionic case. The reason is the larger kaon mass, or equivalently, the smaller kaon 
pair velocity which suppresses emission of such pairs from regions of the hypersurface with 
high flow velocity, i.e., at early times, cf. Fig. 7 (b). (Remember that the emission function 
is large in regions where pair and flow velocity are equal.) This effect leads to considerably 
smaller values for R out and, as already discussed in the context of Figs. 4 and 5, R out /R S ide- 

For larger K\, however, the size of the region in time from which kaon pairs are predomi¬ 
nantly emitted does not decrease very much, as one confirms via the time variance (dashed 
line in Fig. 7 (f)). This at first sight astonishing result is due to the fact that the flow velocity 
remains large up to very late times, cf. Fig. 7 (b), such that kaon pairs with high velocity 
are not only emitted at early times (where the flow velocity is high anyway), but also up 
to late times before the origin cools below T = 0.7 T c . The strong increase in the kaon pair 
velocity with K± leads then to a corresponding increase in R out and R out , cf. Fig. 7 (f).Note 
that due to the larger kaon mass, the difference between dashed and dotted curves is much 
bigger than for the pion case. 
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In Fig. 8 we present the corresponding results for the case AT = 0 which serves to 
explain the behaviour observed in Fig. 6 . In this case, the outward radii increase with K± 
both for kaons and pions. The reason is that the emission functions, besides the shift of the 
maximum to earlier times, do not exhibit a decrease of the emission probability at late times, 
resulting in an overall increase in the duration of particle emission and, correspondingly, in 
the outward radii. From the above it is clear that due to the larger kaon mass, the increase 
of the outward radius with K± is much stronger for kaons than for pions. 

The reason for the increase in the width of the emission function is the behaviour of the 
flow velocity, Fig. 8 (b). Even up to times where the freeze-out reaches the origin, (3 stays 
comparatively large, enhancing the emission probability for particle pairs with high velocity. 
The physical reason for this high transverse velocity is a rarefaction shock wave (cf. Fig. 9 (a) 
in II10II), which travels inwards and expells matter with a constant flow velocity, a fact that 
is confirmed by the plateau observed in Fig. 8 (b) for times between 3 and 6 Rq. Note again 
the excellent agreement between approximate radii, R 0 ut , S ide, and the inverse widths of the 
correlation functions, i? out!Si d e , in Figs. 8 (e,f), although the emission functions are not even 
remotely resembling Gaussian distributions. In fact, for (|16]|T7]) to be good approximations 
of the inverse widths it is sufficient that the correlation functions resemble Gaussians, rather 
than the emission functions themselves. In most (realistic) situations, this happens to be 
the case []. 

We finally mention that the high flow velocity in the cases of interest ensures that the 
negative contributions to the spectra and the correlation functions are negligible, as indicated 
in Subsection 2.3. Quantitatively, they are smaller than 1% (of the total particle number 
crossing the freeze-out hypersurface) for the smallest transverse momenta considered, Ii± ~ 
100 MeV. For larger transverse momenta, they are even smaller. 


4 Summary and Conclusions 


In this work we have investigated the time-delay effect of the QCD transition on two-pion 
and two-kaon correlations in ultrarelativistic heavy-ion collisions in the framework of ideal 
hydrodynamics with transverse cylindrical and boost-invariant longitudinal symmetry. Our 
study is an extension of the work in [Hj, where this time-delay signal, i.e., an enhanced 


ratio R ou t/Rside, was investigated for pion interferometry with fixed average transverse pair 
momentum K± = 300 MeV. Our main interest in the present work was the transverse 
momentum dependence of the proposed signal, and the question whether it can be also seen 
in two-kaon interferometry. Kaon interferometry seems preferable for two reasons. First, 
effects of resonance decays are small (only the shorter-lived K* resonances influence the final 


kaon spectra, at most a (small) 10 %-effect at temperatures corresponding to freeze-out ||T3| ). 
For pion interferometry, however, the decay of long-lived resonances plays an important role 
and virtually increases the (average) source size and lifetime of the system [Q] . This results 
in considerable modifications of the inverse width of the two-particle correlation functions 
of pions |P~5|, |3T|, |55| . Second, neutral kaon correlation functions would not be influenced by 


1 We thank U. Heinz for pointing this out. 
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final state Coulomb interactions at small relative momenta (where the correlation signal is 
large). 

We presented a systematic study of the inverse widths R out , i? S ide, and their ratio for 
two-pion and two-kaon correlation functions at midrapidity, as a function of the initial en¬ 
ergy density eo, the average transverse momentum of the particle pair A'j_, and for different 
equations of state. Studying the emission functions and their variances for selected cases, 
we explained in detail the momentum dependence of the inverse width in out-direction. We 
confirmed that the prolongation of the lifetime of the system due to the QCD transition is 
observable in an enhanced ratio R ou t/R s ide for kaons as well as for pions. The proposed signal 
is maximized for initial energy densities expected to be reached in central Au+Au-collisions 
at the RHIC collider, eo ~ 10 — 20 GeV fm -3 (for realistic [|15j values of the initial ther- 
malization time scale, r 0 ~ 0.5 fm), and for pairs with high average transverse momentum, 
K± 1 GeV. This opens the interesting possibility to observe the time-delay signal of the 
QCD transition with kaon interferometry in the STAR and PHENIX experiments at RHIC. 


Acknowledgements 


We would like to thank M. Gyulassy, T. Hallman, U. Heinz, B. Sclilei, R. Pisarski, and 
W. Zajc for stimulating discussions. D.H.R. thanks in particular U. Heinz for encouraging 
this study, for a critical reading of the manuscript, and for invaluable help in interpreting 
the results. Special thanks go to M. Gyulassy for his continued support and encouragement 
throughout this work. D.H.R. gratefully acknowledges the hospitality of the Institute for 
Theoretical Physics of the University Frankfurt. The work of D.H.R. has been supported 
by the Director, Office of Energy Research, Division of Nuclear Physics of the Office of 
High Energy and Nuclear Physics of the U.S. Department of Energy under contract nos. 
DE-FG-02-93ER-40764 and DE-FG02-96ER40945. 


14 




Appendix 


In this appendix we collect the explicit expressions necessary to evaluate the variances of 
the emission function S(x,K ) appearing in eqs. (p~6|jr7|) . We start with noting that 


(x 2 ) = 

(x 2 ) - (x> 2 , 


(P) = 

<* 2 > - <*> 2 , 

(18) 

xt) = 

(xt) - (x) (t) , 


(?) = 

(?) - (y) 2 , 



while due to the h-function in eq. (O), all space-time variables are taken along the freeze-out 
hyper surface E, given by the 4-vector 


= (Tf K) cosh)), 77(C) cos0, 17(C) si nip, T f ( () sinhi)) , 


(19) 


with t /(C)i r /(C) = tf( (,rj = 0), given by the isotherms in Figs. 7 , 8 (a). Here we adhered 
to the notation of Ref. [[TTJ, i.e., the hypersurface E is parametrized by the three variables 
£, 77, and 0, where ( runs from 0 to 1 along the hypersurface in the (t — r)-plane at z = 0, 
77 from — 00 to 00 along the hypersurface in the (t — z)-plane (77 is identical to the space- 
time rapidity Artanh[z/t]) and 0 from 0 to 2 tt in the transverse (x — 7/)-plane. Due to 
the transverse cylindrical and boost-invariant longitudinal symmetry of the problem, the 
parametrization of the hypersurface in terms of 0 and 77 is trivial. Note that for z = 77 = 0 , 
t = \/t 2 — z 2 = t. 

Using the definition of the average (•} in eq. (0. we find after employing (cf. eq. (B.3) 

u H) 


dE • K = (-H),, cosh 7/ 07 7 K\ coso 0b r/(0 T f (() d(dt)d<f> , 


where we used K° ~ E Kl K = (K±, 0, 0), and (cf. eq. (B.4) of [|T0| ) 

K ■ u — Ek cosh 77 cosh 7] r — K± sinh ? 7 r cos 0 , 


( 20 ) 


( 21 ) 


where rj r = Artanh0 (0 is the radial component of the transverse velocity at z = 0), and 
after integrating over 77 , 0 the final result: 


E 


dr 


Ao(a) + 


(x) = AT 1 d( rjr f |- 

(y) = 0 , 

(t) = A /" -1 dC rfT]{-E K 
(x 2 ) = A f- 1 dC VjTf | -E k ^ A\(a) 

+K ± ^-K 0 {a) 

dC 


dC 

d Tf 

dC 


f K^a) h(b) + K ± ^-K 0 (a) 

dC 


Ad (a) 


m - 


d T f 


h(b) 


,( 22 ) 

(23) 


/o(&)+AA^A0(a)/ 1 (&U,(24) 


m 


h(b) 




(25) 
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(y 2 ) 

(t 2 ) 

(xt) 

J\f 



where a = Ek cosh r) r /T, b = K±smhrj r /T. 


( 26 ) 

(27) 

(28) 
(29) 
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Figure Captions: 

Fig. 1 : (a) The entropy density divided by T 3 (in units of s c /T 3 ), (b) the energy den¬ 
sity divided by T 4 (in units of T c s c /T 4 ) as functions of temperature (in units of T c ), (c) 
the pressure (in units of T c s c ), (d) the square of the velocity of sound as functions of en¬ 
ergy density (in units of T c s c ). The solid lines correspond to AT = 0, the dotted curves 
to AT = 0.1 T c . Quantities for the ideal gas EoS (with d H degrees of freedom) are repre¬ 
sented by dashed lines. The ratio of degrees of freedom in the QGP to those in the hadronic 
phase is dq/du = 37/3. The critical enthalpy density is T c s c ~ 0.75GeV fnT 3 for the case 
dq = 37, du = 3. 

Fig. 2 : -Rout/-R S ide as a function of the initial energy density e 0 for the Bjorken cylinder 
expansion (to = 0.1 Rq) with a first order phase transition (AT = 0). Thin lines correspond 
to pion pairs (m n = 138 MeV), thick lines to kaons (uik = 497 MeV). The pairs are at 
midrapidity and have average transverse momentum K± = 300 MeV. The results are shown 
for freeze-out temperatures Tf = 0.7T C (full lines), Tf = 0.9T c (dotted lines), and Tf = T c 
(dashed lines). 

Fig. 3 : i?out/7?side for kaons as a function of the initial energy density eo at K± = 700 
MeV. The ideal gas case (thin lines) is compared to that of a first order phase transition 
(thick lines, AT = 0). Curves correspond to different freeze-out temperatures as in Fig. 2. 

Fig. 4 : The dependence of R ou t/Rside on the average transverse momentum K± of the 
particle pair, for fixed energy density eo = 18.75 T c s c . Thin lines correspond to pions, thick 
lines to kaons. Different curves correspond to different freeze-out temperatures, as in Fig. 2. 

Fig. 5 : (a) R out /R S ide as a function of eo and K |_ for pions at midrapidity and freeze-out tem¬ 
perature Tf = 0.7 T c , for the case of a first order phase transition (AT = 0, full lines) versus 
the ideal gas case (dotted), (b) as in (a), but for AT = 0.1 T c . (c,d) as in (a,b), but for kaons. 

Fig. 6: (a) Inverse width of the pion correlation function in out-direction R out (in units 
of the initial transverse radius Rq) at midrapidity and for AT = 0, as a function of eo and 
K i_. (b) as in (a), but for R s ic i e . (c,d) as in (a,b), but for kaons. 

Fig. 7: (a) Isotherms in the (t — r)-plane (at z = 0) for the hydrodynamic expansion of a 
Bjorken cylinder with eo = 18.75 T c s c , r 0 = 0.1 Rq, and for a smooth transition AT = 0.1 T c 
in the EoS. Labels correspond to temperatures in units of T c . The freeze-out isotherm 
(Tf = 0.7 T c ) is the thick line, (b) The flow velocity f3 as a function of time along the 
freeze-out isotherm of part (a), (c) The emission function for pions (in arbitrary units), 
integrated over rj and (f> (and smoothed over approximately 100 time steps of the hydrody- 
namical evolution), along the freeze-out isotherm of part (a), for K± = 300 MeV (solid line), 
K l = 700 MeV (dotted line), and K ± = 1000 MeV (long-dashed line), (d) as in (c), but for 
kaons. (e) R^ ut and Tg ide for pions as a function of K± (solid lines), and the approximate 
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radii R^ ut and R 2 ide (dash-dotted lines). The dashed line is the variance in time only, ( t 2 ), 
while for the dotted line, this quantity was multiplied with (i 2 K . For the long-dashed line, 
the x — /-correlation —2 PK(xt) has been added, (f) as in (e), but for kaons. 

Fig. 8 : As in Fig. 7, but for AT = 0. 
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